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Abstract 

The use of conservation laws in nonconservative form for deriving shock jump conditions by 
Schwartz distribution theory leads to ambiguous products of generalized functions. Nonstandard 
analysis is used to define a class of Heaviside functions where the jump from zero to one occurs 
on an infinitesimal interval. These Heaviside functions differ by their microstructure near x = 0, 
i.e., by the nature of the rise within the infinitesimal interval it is shown that the conservation laws 
in nonconservative form can relate the different Heaviside functions used to define jumps in dif- 
ferent flow parameters. There are no mathematical or logical ambiguities in the derivation of the 
jump conditions. An important result is that the microstructure of the Heaviside function of the 
jump in entropy has a positive peak greater than one within the infinitesimal interval where the 
jump occurs. This phenomena is known from more sophisticated studies of the structure of shock 
waves using viscous fluid assumption. However, the present analysis is simpler and more direct. 
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Shock Jump Conditions Using 
Schwartz Distributions 

Assuming that the shock surface is given by 

f(x, t) = 0 and using mass continuity equation in 
conservative form, we have 

+ v • (PH) = ^ + v • (PH) - I V/| [v n Ap - A(p u n )] 5 (/) 

w ,■ '■ v — J 

= 0 =0 =0 Jump Cond. 

d and V are generalized derivatives, A( ) is the 
jump, v = -df/dt/\Vf \ and u = u ■ V/ are shock 

mV mV 

and fluid normal velocities. l~his method fails for 
equations in nonconservative form because 
ambiguous products of distributions appear. 
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Ambiguities of Products of 

Distributions 

Let H(x) be the Heaviside function, then 

H n (x ) = H(x). Taking the derivatives of 
both sides, we get 

nH n ~\x)H'(x ) = H'(x) 

It follows that H n ~ *(x) = l/n. 

Taking n = 2, we get H(x) = 1/2 

OONITriADacTlDiJ! 
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The Source and Solution of 

the Problem 

The jump of the Heaviside function occurs 
at x = 0 which coincides with the support 
of its derivative, the Dirac delta function. 

In fact H(x)S(x) is not defined. This gives 

us a solution to thi_s prob lem 

Let the jump of the Heaviside function occur on 
an infinitesimal interval near x = o and give a 
microstructure to the jump. This rules out 

H n (x) = h(x) and we have a new class of 
Heaviside functions each different from others. 
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A New Class of Heaviside Functions 



The jump of the Heaviside functions occurs on a true 

infinitesimal interval in the halo of x = 0. We, thus, have 
a class of Heaviside functions all basically behaving as 
the classical Heaviside function and all satisfying 

Yi 

H ( )*m x ) . We claim that the conservation laws in 
nonconservative form give the relations between the 
various Heaviside functions associated with the jump. 
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Some Properties of the New 
Class of Heaviside Functions 

i) Two Heaviside functions are equal 
if they have the same microstructure 
and infinitesimal interval of the jump. 
The microstructure is assumed to be 
as smooth as needed in a given 
problem. 

ii) The derivative of a Heaviside 
function behaves as the Dirac delta 
function. The derivative takes 

hyperreal infinite values on the 
infinitesimal interval of the jump. 
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Some Properties of the New Class 
of Heaviside Functions (Cont’d) 

iii) The usual algebraic manipulations 

of smooth functions can be applied 
to the Heaviside functions in the new 
class. 

iv) Product of a Heaviside function 
h(x) and its derivative h\ X ) is always 
defined. 

v) Cancellat on law holds: ew and pm 

smooth, then h^ u \x)q(x) = h^ , 1 \x)p(x) im plies 
that q(x) = p(x)m 
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Shock Jump Conditions- 
Nonconservative Form 


Assumptions: 

i) Different Heaviside functions 
H(^), K(%), l(£) and r(£) describe the 
jumps of specific volume v, velocity u, 
pressure p and entropy s, respectively. 

ii) The jump for all flow parameters 
occur on the same infinitesimal 

interval (o, e) . 
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Shock Jump Conditions- 
Nonconservative Form (Cont’d) 

Let \ = x-ct and A( ) = ( ) 2 - ( the jump in 

a fluid parameter. Define 


v = Vj + Avif(^) Specific volume 


u = u { + AuK(^) Fluid velocity 

Shock 

p = p l + ApL(b > ) Pressure 

► c 

(i) 

S = S l +AST(t ) ) Entropy 

(2) 

c is the shock speed 
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Shock Jump Conditions: 
Nonconservative Form (Cont’d) 


Mass conservation law v, - vu + uv = 0 

t X X 


gives — c)AvH' — v^AuK'- AvAuHK' + AvAuH'K = 0 


resulting in the ODE 


X U] - c 

K + 


V 


A u 


\dH _ v _\_ 
JdK Av 


V, / H, 

The solution is H{%) = - — + b\K + — 


n, - c 


\ 


Av 


Am J 


, b const. 


From H( 0) = K( 0) = 0, we get b = v.Au/[Av(u. -c)] 


and from //(e) = K(e) = 1, we get *> = l . It 

follows that H{t,) = K( and c = (u^ Av - Au)/Au . 
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Shock Jump Conditions- 
Nonconservative Form (Cont’d) 


Momentum equation m. + mm+vz? = 0 gives 

( A A 


/• M 1 — C 

K + 


V 


Am 


2 / v i \ 

(Am) K' + ^ + — jAvApL' = 0. We have 


used //((;) = K(%) here. This gives the ODE 


dL (Am) _ Q. We used (u,-c)/Au = v,/Av 
z/A' AvAp 

from fo = l. The solution of this ODE is 

2 

Lit,) = + const.. From L(0) = ^(0) = 0, 

we get d = o. From L(e) = /T(e) = 1 , we get 

(Am) 2 / AvA/? = -l. It follows that L(^) = A^). 
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Shock Jump Conditions- 
Nonconservative Form (Cont’d) 


The constitutive relation S = c v \og(pv y ) 

gives S = S 1 +AST{fy = c log[(l + ApH)(l + AvH) y ] 


= c v lo g [ Pl viyc v lo g 


( A v Y Av Y 

1 + ^-H\ 1 + — H 

V Py A v 


1 


Using s y = 


c v iog(p jv|) , we get 


Nonlinear ! 
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The Behavior of T(x) 


If we assume that the Heaviside function 
H(x ) is monotonic, then T(x ) will have a 
maximum above 1 within the shock layer 
as predicted by more complex analysis. 
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Summary of the Results 


We have shown that 


H(Z>) = K&) = L(^) 




( 

1 + 

v 




Shock speed: c - (u^Av - v^Au)/Au 


The relation ( Au) 2 /AvAp = -l gives the 3 randtrs 
relation (m, -c)(m 2 -c) = Ap/Ap. All other known 

shock jump conditions are also obtained. 
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Concluding Remarks 


— We defined a new class of Heaviside 
functions with smooth jump on an 

infinitesimal interval in the halo of x = 0 and 
satisfying H n (x)*H(x). 

— Using this class of Heaviside functions, we 
obtain the shock speed and jump conditions 
from conservation laws in nonconservative 
form. 

— The relative shape of the microstructure of 
the jump of different fluid parameters can be 
obtained by rigorous algebraic manipulations. 

— The unusual behavior of entropy within a 
shock can be obtained by using the new class 
of Heaviside functions. 
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